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The existence of macroscopic spin currents in the ground state of superconductors is predicted 
within the theory of hole superconductivity. Here it is shown that the electromagnetic Darwin 
interaction is attractive for spin currents and repulsive for charge currents. It is also shown that 
the mere existence of spin currents implies that some electrons are moving at relativistic speeds in 
macroscopic superconductors, which in turn implies that the Darwin interaction plays a fundamental 
role in stabilizing the superconducting state. 
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The theory of hole superconductivity0| predicts that 
superconductors expel negative charge from the inte- 
rior towards the surface 0- As a consequence an excess 
of negative charge is predicted to exist within a Lon- 
don penetration depth of the surface and excess positive 
charge in the interior. A description of superconductors 
as 'giant atoms' results|3j, and an alternative electrody- 
namic theory Hwhich is different from the conventional 
London theorvQ and leads to experimentally testable 
consequences |g. 

As a consequence it is predicted that macroscopic spin 
currents should exist in the superconducting state 0, 
in the absence of applied fields, as shown schematically 
in Figure 1. Classically this can be understood from the 
requirement that electrons near the surface should move 
in macroscopic orbits of large angular momentum so that 
the centrifugal force prevents them from 'falling in' to- 
wards the interior to neutralize the positive charge. Be- 
cause in the absence of applied magnetic field there is no 
charge current in the superconductor, electrons of oppo- 
site spin will orbit in opposite direction. The preferred 
direction is determined by spin-orbit coupling, i.e. the 
electron magnetic moment points predominantly parallel 
to its orbital angular momentum. Thus in the absence 
of applied fields Cooper pairs (fc t, —k I) carry no charge 
current but carry spin current. 

Superconductivity is of course a quantum phe- 
nomenon. However the purpose of this paper is to show 
that classical arguments can shed important light into 
the fundamental physics of superconductors. In the sense 
of Bohr's correspondence principle we argue that the 
macroscopic quantum manifestations of superconductiv- 
ity should also be understandable from a classical point 
of view. Here we focus on the Darwin interaction, which 
is the lowest order relativistic correction to the electrody- 
namic Lagrangian of interacting charged particles 

That the Darwin interaction plays a fundamental role 
in superconductivity was first proposed in prescient work 
by H. Essen in 1995 llOl. and advocated by him in sev- 
eral papers thereafter |llL IT^ . Essen presents a variety 
of microscopic and macroscopic arguments in support of 
his assertion that the Darwin interaction leads to an at- 




FIG. 1: Schematic picture of Cooper pairs in a spherical su- 
perconductor giving rise to spin currents. The arrow parallel 
to the orbit indicates the direction of motion and the arrow 
perpendicular to the orbit indicates the direction of the elec- 
tron magnetic moment to give minimum spin-orbit energy. 
Excess negative charge exists within a London penetration 
depth of the surface, and excess positive charge in the inte- 



tractive interaction between electrons and consequently 
to Cooper pairing and to superconductivity. However, 
in Essen's treatment the Darwin interaction is found to 
be attractive for parallel currents, hence he predicts the 
existence of charge currents in the ground state of su- 
perconductors, similar to old theories of Heisenbergjl^j 
and of Born and Cheng[T^. which is however proven 
impossible by 'Bloch's theorem' [l5l|. According to Es- 
sen's treatment the interaction would be repulsive for 
spin currents, and in contrast to our work there is no 
charge inhomogeneity in his description of superconduc- 
tors. Furthermore, in Essen's treatment electrons move 
at non-relativistic speeds and hence the magnitude of 
the Darwin energy is small, and certainly cannot account 
for high Tc superconductivity ^3 ■ Instead we argue here 
that the interaction is in fact attractive for spin currents 
and repulsive for charge currents, and in addition that 
electrons in spin currents move at relativistic speeds. In 
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our work, no charge currents exist in superconductors in 
the absence of apphed fields, and 'Bloch's theorem' does 
not preclude spin currents. 

The possible role of the Darwin interaction in super- 
conductivity was also considered in the work of Capelle, 
Gross and coworkers |16|. However in that work it was 
not found that the Darwin term plays an essential role in 
the attractive interaction leading to pairing, rather rel- 
ativistic corrections were considered 'independent of the 
interaction itself. The role of the Darwin interaction 
in superconductivity was also considered in the work of 
Lurie and CremerjlTj. however its possible role in the 
pairing interaction was not discussed there either. 

There has been confusion in the literature about the 
sign of the Darwin interaction. In applications to plasma 
physics, early work used an approximate form that 
predicts an attractive interaction for parallel (charge) 
currents[3|- This was later corrected"]^ and it is now 
agreed by several workers that the interaction is in fact 
attractive for antiparallel currents[ll I20ll2lll2^ , while 
others continue to argue that the opposite is true[23. 

The Darwin electromagnetic interaction Lagrangian 
between two charges qi and q2 includes the lowest order 
relativistic correction and is given byQ 
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[~l + ^Wi ■ V2 + {vi ■ f)iv2 ■ r)]] (1) 



where vi , V2 are the velocities and f is the unit vector con- 
necting the positions of both charges. The Hamiltonian 
is obtained from the Lagrangian through the relation[2^ 
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with 
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hence for a time-independent situation where the Hamil- 
tonian gives the constant energy of the system, the inter- 
action energy is 

E.nt = ^ [1 + ^ [«1 • + (t/l • f) {V2 ■ f)]] (3) 

Note that in the usual situation where L = T ~ V ^ H — 
T+V, with T and V the kinetic and potential energy, the 
sign of the second term in Eq. (1) would have changed 
in going from Eq. (1) to Eq. (3). This is not the case 
here because the potential is velocity-dependent. 

Inspection of Eq. (3) reveals that the interaction en- 
ergy between two charges is lowered when their velocities 
have opposite sign, and that the magnitude of the sec- 
ond term becomes important compared to the ordinary 
Coulomb repulsion for speeds approaching the speed of 
light. Both of these facts are key to the role of the Darwin 
interaction in superconductivity within our theory. 



The confusion about the sign of the Darwin interaction 
arises because expression (3) appears to contradict the 
well-known fact that in magnetostatics parallel charge 
currents attract. Hence it is argued that Eq. (3) is 
'misleading' 23J and instead Eq. (1) should be inter- 
preted, using L = T — V^, as representing a potential 
energy that is lower for parallel currents. In fact, as dis- 
cussed by Schwingerj2^ the correct expression for the 
interaction energy in terms of charge and current density 
for a set of charges qi with velocities Vi is 
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_L j d\d\'^^p{r,t)^^p{r^,t)\r-7^\ (4) 



with 



and using the continuity equation 



dp 
dt 
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Eq. (3) results for each pair of particles. The second term 
in Eq. (4) also appears to contradict the expectation that 
the energy should be lower for parallel currents, however 
when combining the second and third terms in Eq. (4) 
Schwinger shows that they reduce to 
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predicting lower energy for parallel currents. However in 
obtaining Eq. (7) from Eq. (6) a boundary term, that is 
indeed irrelevant for long parallel wires, is dropped. For 
our case, Eq. (7) is not applicable. 

The following physical argument may help clarify the 
situation. The interaction between parallel wires with 
currents of the same sign is attractive because the mag- 
netic field generated by the wires has opposite sign in the 
region between the wires. The magnetic field is strongest 
close to the wires, and because it contributes to the mag- 
netic energy through 
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it will lower the energy to bring the wires closer together 
to cancel the magnetic fields in the region close to the 
surface of each wire facing the other wire. For the same 
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reason, the force between parallel wires with antiparal- 
lel currents is repulsive. However, consider two identical 
current loops made of wires of finite cross section so that 
the planes of the loops are parallel and their centers lie on 
a common axis perpendicular to the planes of the loops, 
and allow to change their relative distance d including the 
regime where the loops overlap. For parallel currents, the 
magnetic energy for d = is twice as large as for d = oo 
(four times as large as the energy of a single loop, be- 
cause the current is twice as large), while for antiparallel 
currents the magnetic energy when d = is zero because 
the currents have cancelled out. Consequently the in- 
teraction in the case of antiparallel currents has to be 
attractive in some regime (since the energy gets lowered 
from finite to zero as d goes from oo to 0) and for the case 
of parallel currents it has to be repulsive in some regime, 
to increase the energy. Clearly that regime is when the 
loops are close and the currents overlap. We conclude 
that for charged particles circulating in the same region 
of space indeed the energy will be lowered when the par- 
ticles circulate in opposite directions and no net charge 
current exists so that no magnetic field is generated. 

For two electrons of opposite spin orbiting in opposite 
direction in a classical orbit of radius R the energy Eq. 
(3) is 

r 4 

where 9 is the angle between the position vectors of both 
charges which changes with time according to 9 = 2u)t, 
uj = V / R, and r is the instantaneous distance between 
the charges, r — \/2R^/l — cos9. Note that when 9 
the second term in Eq. (9) exactly cancels the first 
term, the ordinary Coulomb repulsion, as v c. For 
any angle < 9 < 109.5° the velocity-dependent term 
is negative and reduces the Coulomb repulsion. Small 
values of 9 correspond of course to small r which is when 
the ordinary Coulomb repulsion is largest. 

In the model of hole superconductivity a uniform posi- 
tive charge distribution exists in the interior of supercon- 
ductors giving rise to a radial electric field E = cr, with 
c of order lO^V/ cm'^W\. The angular velocity of electrons 
is then uj = (ec/mg)^/^ with rrie the electron mass, and 
the second term in Eq. (9) will become very important 
for macroscopic samples. 

In a quantum- mechanical treatment, assuming a spher- 
ical geometry the Cooper pair can be denoted by 

with n, I, m radial, orbital and azimuthal quantum num- 
bers. It represents electrons orbiting in opposite direc- 
tions as in the classical example discussed above, hence 
we conclude that the energy of a Cooper pair is lowered 
by the velocity-dependent Darwin interaction. There is 
yet another relativistic lowering of the energy of Cooper 
pairs, from the spin-orbit interaction with the internal 



electric field E 

Es.o. = TT^S ■ {v X E) (11) 

with S the electron spin. Even though this term goes as 
(w/c) compared to the Darwin term's (v/c)^ dependence, 
the prefactor is small and at relativistic speeds the energy 
lowering from this term is negligible compared to the 
Darwin term. 

The expression Eq. (9) diverges as ^ 0, hence cannot 
be used to estimate the energy. As a simple example, 
consider two spherical shells of radius R and charge q 
rotating in opposite direction with angular velocity co, 
corresponding to the electrons with spin 'up' and 'down' 
respectively. In a stationary situation the third term in 
Eq. (4) vanishes and we can calculate the Darwin energy 
from just the second term in Eq. (4) 

Eoar^^n = ^J d^rd^r'^J^^ (12a) 

with 

and cr = ±1, and obtain 

q'co^R_ NWR 

Ji^Darwin - g ^2 ~ g^2 V"^"^^ 

for q = Ne, with N the number of Cooper pairs. The 
kinetic energy is (neglecting relativistic corrections) 

Ek^n = ^Nm^uj^R^ (14) 

so that for large N the Darwin energy lowering is greater 
than the kinetic energy associated with the spin current. 
The relevant parameter (ratio of Eqs. (13) and (14)) is 
Ne^ /nieC^R as already found by Essen 12). For relativis- 
tic speeds the energy lowering Eq. (13) is of the same 
order of magnitude as the Coulomb energy cost of the 
inhomogeneous charge distributional. 

Note also that the rigid rotation assumed in Eq. (12b) 
is far from optimal to benefit from the Darwin energy 
lowering, since only the electrons in the equatorial plane 
are moving at the maximum speed v = uR. Other pat- 
terns of motion where electrons move along various great 
circles of different orientations as depicted in Fig. (1) 
will yield larger energy lowering than Eq. (13). In ref. 
|2^ , a numerical simulation of electrons on the surface of 
a sphere was performed and it was found that electrons 
spontaneously develop complicated motion patterns that 
are indeed characterized by neighboring electrons having 
velocities in opposite directions. 

In the theory of hole superconductivity, a negative 
charge density p_ exists within a London penetration 
depth of the surface, and a positive charge density po in 
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the interior. The charge density obeys the differential 
equation 



p{r) = po + AiVV(f) 



(15a) 



with Al the London penetration depth, and the electric 
field is determined by the equation 
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where Eq is the electric field originating in a uniform 
charge density po- We consider a spherical geometry for 
simplicity, where Eq is given by 



(15c) 



and q — AttR^pq/S is the charge expelled from the inte- 
rior of the sphere of radius R to the surface shell. The 
Darwin interaction term becomes important compared 
to the ordinary Coulomb repulsion when the speed ap- 
proaches the speed of light. We now show that in fact the 
predicted spin currents involve such speeds for samples 
of dimensions much larger than the London penetration 
depth. 

Consider a sphere of radius R >> Xl, and electrons 
associated with the spin current near the surface rotating 
with azimuthal velocity vq, in opposite directions for spin 
1 and |. Mechanical equilibrium requires that 
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for electrons at radius r. We will use nonrelativistic ex- 
pressions for simplicity since it is sufficient to illustrate 
the point. When a magnetic field B is applied, the speed 
of electrons within a London penetration depth of the 
surface is modified to screen the magnetic field, denote 
the change in speed by v^. To maintain mechanical equi- 
librium to first order it is seen from variation of Eq. (16) 
that the condition 
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is required. Now from integrating the London equation 



V X J= - 



B 



we can write the charge current J as 
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and defining the current density as 
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we conclude that the charge density p that contributes 
to the current is given by 
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where r i? to within a London penetration depth cor- 
rection. Using for the London penetration depth 
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with Us the density of electrons in the condensate, Eq. 
(20) yields 



p = 2ens — 
r 



(22) 



which differs in a fundamental way from the usual as- 
sumption in London theory. Namely, it says that only a 
small fraction ~ (Xl/t) of all electrons contribute to the 
transport current, with a speed that is larger by a fac- 
tor r/Ai from the speed assumed in conventional London 
theory. 

The maximum value that the electric field attains is 
Emax q/R^ ~ AttXlP- within distance Xl of the 
surface^. Here, q is the expelled charge and p_ is the 
excess negative charge density near the surface. This 
electric field determines the maximum speed of the spin 
current vq (Eq. (16) for r ~ R). Under the reasonable 
assumption that the transport charge density p given by 
eq. (20) is of the same order of magnitude as the ex- 
pelled charge density one is led from Eqs. (16) and 
(20) to the startling conclusion tha t fn ~ c in samples 
of dimensions much larger than Ar |26j|. This then im- 
plies that the Darwin interaction is all-important in the 
energy balance. 

In addition, the condition p = p^ leads through Eq. 
(20) to the remarkable conclusion that 
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in other words, that the potential energy of electrons 
within Al of the surface where the electric field is max- 
imum is right at the threshold where pair production 
becomes energetically possible! . Possible experi- 

mental consequences of this finding will be discussed else- 
where. 

Finally we may ask, how many Cooper pairs partic- 
ipate in the charge and spin currents near the surface? 
Denoting by N the number, we have Ne — 2ttR^XlP, 
and using Eq. (20) for p 
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with Ve the classical electron radius. This is precisely the 
condition derived by Essen [T^ for the Darwin interac- 
tion to be important. Furthermore, as already showed 
by Essen 12], it is the condition required to explain the 
magnitude of magnetic fields in rotating superconduc- 
tors (London field) under the assumptions that only a 
fraction n^Ai/i? of electron density (Eq. (22)) near the 
surface contribute to the current and that these electrons 
are completely unaffected by the body' s rotation (these 
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assumptions were not made explicit in Essen's work). In 
contrast, in the conventional London treatment @ all the 
Us electron density near the surface contributes to the 
London field through a velocity 'lag' Av ^ ojXl for each 
carrier, with lu the angular velocity of body rotation. 

In summary, we have given a series of arguments that 
lead to the conclusion that relativistic effects play a fun- 
damental role in superconductors of dimensions larger 
than the London penetration depth within the theory of 
hole superconductivity. This was already foreshadowed 
by the prediction^ that macroscopic electric fields of 
magnitude of order lO^V/ cm exist in the interior of su- 
perconductors, which implies relativistic effects in sam- 
ples of dimension of order cm (recall that the electron 
rest energy is ~ 0.5 x lO^eV). Separately, the mere 
existence of charge expulsion and spin currents and the 



condition of mechanical equilibrium was shown here to 
imply that electrons near the surface move at relativistic 
speeds, which in turn implies that the Darwin electro- 
magnetic interaction has to play an important role. We 
have argued that the Darwin interaction is attractive for 
spin currents and repulsive for charge currents (in con- 
tradicion to previous work), and hence conclude that it 
plays a fundamental role in the stabilization of the su- 
perconducting state in the theory of hole superconduc- 
tivity where spin currents are predicted to exist. The fact 
that the Darwin energy lowering is a kinetic energy low- 
ering mechanism was also foreshadowed by the lattice 
formulation of the model of hole superconductivity [23. 
Further development of these concepts will be given in 
future work. 
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